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1. Introduction 

We survey results of recent activity aimed at studying controllability and 
accessibility properties of Navier-Stokes systems (NS systems) controlled by 
low-dimensional (degenerate) forcing. This choice of control is characteristic 
feature of our problem setting. The corresponding equations are 

(1) du/dt + V u « + gradp = uAu + F(t, x), 

(2) divu = 0. 

The words "degenerate forcing" mean that F(t,x) can be represented as: 

F(t,x) = v k (t)F k (x), /C 1 is finite. 
fee/c 1 

The word "controlled" means that the functions Vk(t), t £ [0, T] entering the 
forcing can be chosen freely among measurable essentially bounded functions. 
In fact any functional space, which is dense in Li[0, T] would fit. 

The domains treated here include 2-dimensional (compact) Riemannian 
manifolds M homeomorphic to either sphere, torus or disc. The latter case 
includes a rectangle and 2D simply connected domain M with analytic 
boundary dM. We impose so called Lions boundary condition, whenever 
boundary is nonempty. 

Our approach stems from geometric control theory, which is essentially 
based on differential geometry and Lie theory; geometric control approach 
proved its effectiveness in studying controlled dynamics in finite dimensions. 
We would like to report on some ideas of how such methods can be extended 
onto the area of infinite-dimensional dynamics and of controlled PDE. Ex- 
tensions of the geometric control theory onto infinite-dimensional case are 
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almost unknown. Classical Lie techniques are not well adapted for infinite- 
dimensional case, and several analytic problems are encountered. 

In this contribution we concentrate almost exclusively on geometric and 
Lie algebraic ideas of the accomplished work. For details on analytic part 
we refer interested readers to the publications [3 EH EH ESI l28l 129] . 

Applications of geometric theory to the study of controllability of finite- 
dimensional systems is well established subject, although many problems 
still remain unsolved. Starting point of the activity aimed at controlling NS 
systems by degenerate forcing was study (p~H IH [6j El]) of accessibility and 
controllability of their finite-dimensional Galerkin approximations on T 2 and 
T 2 (periodic boundary conditions). This question being settled, one should 
note, that controllability of finite- dimensional Galerkin approximations of 
NS systems on many other domains remains an open problem; answers for 
generic analytic 2D domains follow from results of Section [TUl 

Study of infinite-dimensional case started in [HEIE], wnere we dealt with 
2D NS/Euler system on 2D torus T 2 . In those publications notions of solid 
controllability in projections and of approximate controllability have been 
introduced and sufficient criteria for them have been established. 

To arrive to such criteria the technique of so-called Lie extensions in in- 
finite dimensions has been suggested. In the context of our problem this 
technique can be loosely interpreted as designing the propagation to higher 
modes of the energy pumped by controlled forcing into lower modes. 

The control functions involved are fast-oscillating and analytic part of 
the study consists of establishing continuity properties of solutions of NS 
systems with respect to so called relaxation metric of forcing. The latter 
metric is weaker than the classical metrics and is adapted for dealing with 
fast oscillating functions. 

An extension of the above mentioned techniques techniques onto the case 
of NS system, subject to Lions boundary conditions on a rectangle, has 
been accomplished by S.Rodrigues ([24]). In the course of this study both 
geometric and analytic part needed to be adjusted: Lie extensions turn more 
intricate and continuity properties need to be reproved. These results are 
surveyed in Section [U 

An additional analytic effort (perturbation techniques) is needed for es- 
tablishing controllability criteria for NS system on smoothened rectangular; 
this is described in Section [H 

A new approach is suggested for establishing controllability on a generic 
2D domain (Section [THJ). 

Finally study of Lie algebraic properties of spherical harmonics results in 
controllability criterion for NS/Euler system on 2D sphere (Section [TT1). 

The results appearing in Sections l9l[TT1 have not been previously published. 

An interesting extension of the above described methods onto the case of 
NS system on 3D torus has been accomplished by A.Shirikyan in [28j [29] . 
The geometric part of his study coincides essentially with the one of [6] and 
of [27], but additional analytic difficulties in the 3D case are numerous. We 
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do not survey these results for the sake of remaining more geometric in spirit 
and avoiding new notation. Interested readers should consult [28, 29]. 

There was an extensive study of controllability of the Navier-Stokes and 
Euler equations in particular by means of boundary control. There are var- 
ious results on exact local controllability of 2D and 3D Navier-Stokes equa- 
tions obtained by A.Fursikov, O.Imanuilov, global exact controllability for 
2D Euler equation obtained by J.-M. Coron, global exact controllability for 
2D Navier-Stokes equation by A.V. Fursikov and J.-M. Coron. The read- 
ers may turn to the book [15] and to the surveys [16] and [H] for further 
references. 

The authors are grateful to S.Rodrigues for useful comments and help 
during preparation of this contribution. 

2. 2d ns/euler system controlled by degenerate forcing. 
Definitions and problem setting 

2.1. NS/Euler system on 2D Riemannian manifold. Representation 
of NS/Euler system in the form ([TJ) - (J2J) requires interpretation, whenever 
one considers the system on 2D domain M with arbitrary Riemannian met- 
ric. There is a general way of representing NS/Euler systems on any re- 
dimensional Riemannian manifold (see e.g. |10|). but we prefer to remain 
in 2 dimensions and to advance with some elementary vector analysis in the 
2D Riemannian case. 

We consider smooth (or analytic) 2-dimensional Riemannian manifold M 
(with or without boundary), endowed with Riemannian metric (•, •) and 
with area 2-form a. All functions, vector fields, forms, we deal with, will be 
smooth. 

Any vector field y on M can be paired with two differential 1-forms defined 

as: 

(3) y^y*: {vt,Z) = <T(v,0, 

for each vector field £. Obviously = a(y,y) = 0. 

Note that for any 1-form A there holds: 

(4) A Ay' = (A, y)a. 

To prove {I} it suffices to compare the values of 2-forms A Ay^ and (A, y)a on 
any pair of linearly independent vectors. Evidently (@J) is valid if y (and y$) 
vanishes. If y ^ 0, we take a pair y, z, which is linearly independent. Then 

(\Ay i )(y,z) = 

Now we define the vorticity curl and the divergence div of a vector field 
via the differentials dy ,dy$. These latter are 2-forms; we put 

dy = (curly)cr, dy^ = (divy)a, 



(X,z) 



(y\y) 
(yKz) 
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or by abuse of notation: 

(5) (curly) = dy b /a, (divy) = dy^/a. 

The gradient grad(/j of a function <p is the vector field paired with dtp 
metrically: (grad^)^ = dip. 

As in the Euclidean case the vorticity of gradient vector field vanishes: 

curl(grad9j) = d(gr&dipf / a = d(dip)/a = 0. 

While in 3D case operator curl transforms vector fields into vector fields, 
in the 2D case it transforms vector fields into scalar functions (actually in 
a component of vector field directed along additional third dimension). We 
will define the vorticity operator curl on functions. The result of its action 
on a function (ft is a vector field curie/), which satisfies the relation: 

(A, curlier = (d(ft A A), 

for each 1-form A. By virtue of (J4j) and due to the nondegeneracy of paring 
y i ^ y we conclude: 

(6) (curlc/>) tt = -d(ft. 

As in the Euclidean case the divergence of vorticity of a function vanishes: 

(7) div(curl0) = d{curl(ftf/a = -d(d(ft)/a = 0. 

Coming back to the equation J2]) we note that the condition divu = can 
be written down as 

(8) dJ = 0. 

If M is simply connected we conclude that it" must be a differential: it" = 
—dift, where tft is so called stream function. By virtue of j6j) 

curl?/; = u. 

For non simply connected domains we impose a condition which guarantees 
the exactness; in the next subsection we comment on it. 

Given the symplectic structure on M defined by a and (■, •) we see that u 
is Hamiltonian vector field corresponding to the Hamiltonian — tft: u = — ift . 

The nonlinear term S7 u u in the right-hand side of ([TJ) corresponds to co- 
variant derivative of the Riemannian (metric torsion-free) connection on M. 

Finally we define the Laplace-Beltrami operator A as 

A = curl 2 . 

In Hodge theory (see [10J ) this operator transforms p-forms into p-forms; in 
our notation A transforms vector fields into vector fields and functions into 
functions. 
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2.2. Helmholtz form of 2D NS sytem. To arrive to the Helmholtz form 
of the NS system ([I])-([2j), we apply the operator curl to both parts of (JTJ) . 
As a result we get for the vorticity curlu = w the equation 

(9) dw/dt + curl(V M u) = uAw + f(t, x), 

where f{t,x) = curLF(i, x). 

One should note that the vorticity of gradp vanishes and that the operator 
curl commutes with A = curl 2 . 

To calculate curl(V u ti) according to the formula Jl])we first compute the 
1-form (V u u)^, adapting the argument of pJJl §IV.1.D]. 

Let y be a vector field which commutes with u: the Lie-Poisson bracket 
[it, y] = 0. Then 

(10) ((V M u) b , y) = (V u u, y) = L u {u, y) - (u, V u y). 

(Here and below L u denotes Lie derivative. Note that for the covariant 
derivative of metric connection there holds: L u (u, y) = (V u it, y) + (u, V u y).) 
Since the connection is torsion-free and [u, y] = 0, then V u y — V y u = 0, and 
the right-hand side of (flOl) can be represented as 

L u (v*,y) - (u, V y u) = L u (u\y) - ^{d(u, u),y). 

Besides L u (u b ,y) = (L u u^,y), as long as L u y = [u,y] = 0, and we obtain 

((X7 u uf,y) = {L u u\y) - -(d(u,u),y). 

As far as one can find vector field y, which commutes with u and has any 
prescribed value at a given point, we conclude: (V^u) 1, = L u u b — hd(u,u). 
Using the definition of curl j5]) we get 

curl(V u u) = d({S/ u uf) / a = dL u vt ' /a = L u du 'jo = L u (wa)/a = L u w. 

Hence curl(V u n) = L u w. 

For u being Hamiltonian vector field with the Hamiltonian — ip there holds 
V u w = —{tp,w}, where {•, •} is Poisson bracket of functions. 

The Helmholtz form of the Navier-Stokes equation (cf. [10]) reads 

dw r ; -i a rt \ 

— - {ip,wj - uAw = f{t,X). 
Note that w = curlu = cwc\ 2 ip = Aip. 

The Lions condition written in terms of the vorticity w and the stream 
function tp reads: 

( U ) ^\dM = W \ d M = °- 

If the boundary dM of M is smooth, then the Hamiltonian vector field 
u = — ifi is tangent to dM. 

Given the vorticity w and the boundary conditions (fTTj) . one can recover 
in a unique way the velocity field u, which corresponds to an exact 1-form 
uK The corresponding formula is u = curl^, where tp is the unique solution 
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of the Dirichlet problem Aip = w under boundary condition (jlll) . Indeed 
such u is divergence free and its vorticity equals w by the definition of A. 
The NS system can be written as 

(12) ^ - {A- 1 ™, «;} - uAw = f(t,x). 

This last equation looks universal: in fact its dependence on the domain 
is encoded in the properties of the Laplacian A on this domain. It is well 
explained in [8j [10] that the Euler equation for fluid motion is infinite- 
dimensional analogy of Euler equation for rotation of (multidimensional) 
rigid body, and the Laplacian in (TT2l) plays role of tensor of inertia of rotat- 
ing rigid body. 

2.2.1. Stream function on flat torus. Let us consider the flat torus T 2 , en- 
dowed with standard Riemannian metric and with area form a, both inher- 
ited from covering of T 2 by the Euclidean plane. Let ipi, if2 be the 'euclidean' 
coordinates on T 2 . We proceed in the space of velocities u with vanishing 
space average: J T2 uda = (due to flatness we may think that all the veloc- 
ities belong to the same linear space). 

To establish exactness of the closed 1-form vr (involved in (JHJ)) it suffices 
to prove that its integral along a generator of a torus vanishes. By Stokes 
theorem the integrals of the closed form it* along any two homologous paths 
have the same value. 

Taking u = (ui, 112) we get u* = —u%dipx + U\dif2, Integrating w along a 
loop T : (pi = a, we obtain the value of the integral f T v$ = J* Q 27r u\d<^i = 
c(a), which by the aforesaid is constant: c(a) = c. Integrating it with 
respect to cpi we conclude 2irc = J Q 27r u\(kp%dip\ = f T2 U\da = 0. Hence 
/ r u' = c = 0. The same holds for the loops V : <p% = const. 

2.3. Controllability: definitions. In what follows we reason in terms of 
so called modes which are the eigenfunctions <j> k {x) of the Laplace-Beltrami 
operator A defined in the space of vorticities w: A(j) k (x) = Xk4> k (x). 

Representing w,f in (fl~2l) as a series w(t,x) = Y^k { ?fc(^) l ^ fc (^), f(t,%) = 
Ylk v ^(^)'i )k ( x ) with respect to the basis of eigenfunctions one can write the 
NS system as an infinite system of ODE on the coefficients qk(t). Assume 

W(x)^(x)} = Y,C l ^ k {x). 

k 

Then the equation (fT2l) can be written down 'in coordinate form' as 

(13) C]^\~ l qiqj - v\ k q k = v k {t). 

Typically we will consider the controlled forcing, which is applied to few 
modes 4> k (x), k G K\ K, 1 - finite. Then in the system Jl3j) the controls enter 
only the equations indexed by k € JC , while for k JC \ v k = 0. 
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Introduce another finite set K° of observed modes. We will always assume 
K° D JC 1 . We identify the space of observed modes with ~R N and denote by 
n° the operator of projection of solutions onto the space of observed modes 
span{0fc| k € 1C°}. The coordinates corresponding to the observed modes 
are reunited in observed component q°. 

Galerkin /C°-approximation of the 2D NS/Euler system is the equation 
(ODE) for q°(t), obtained by projecting the 2D NS system onto the space 
of observed modes and putting all the components qk(t), k £ JC° zero. The 
resulting equation is 

(14) -jL - WiA^q^q"} - vAq° = f(t,x). 

As far as K° D K , i.e. controlled forcing / only affects part of observed 
modes, then Wf(t,x) = f(t,x). 

In coordinate form passing to Galerkin approximation means omitting the 
equations (fT3l) for variables qk with k £ JC° and taking these qk equal in 
the resting equations. 

We say that a control f(t,x) steers the system (Qjj) (or [14\)) from (p to 
tp in time T, if for equation (fT2j) forced by / the solution with the initial 
condition (p at t = takes 'value' (p at t = T. 

The first notion of controllability under study is controllability of Galerkin 
approximation. 

Definition 2.1 (controllability of Galerkin approximation). Galerkin K°- 
approximation of 2D NS/Euler systems is time-T globally controllable if for 
any two points q,q in , there exists a control which steers in time T this 
Galerkin approximation from q to q. □ 

This is purely finite-dimensional notion. The next notion regards finite- 
dimensional component of solutions, but takes into account complete infinite- 
dimensional dynamics. Let us introduce some terminology. 

Definition 2.2 (attainable sets of NS systems). Attainable set of the 
NS/Euler system 111 S\) is the set of points in H 2 (M) attained from <p by 
means of essentially bounded measurable controls in any positive time. For 
each T > time-T (time-< T) attainable set (A^ T ) of the NS/Euler 
system is the set of points attained from (p by means of essentially bounded 
measurable controls in time T (in time <T). Attainable set Ag> = Ut-^- '— ' 

Definition 2.3. The NS/Euler system is time-T globally controllable in pro- 
jection onto C if for each (p the image Il £ f^^l coincides with C 

Definition 2.4. The NS/Euler system is time-T L2- approximately control- 
lable if A 1 - is L2-dense in H 2 . □ 

Let us introduce the notion of accessibility in projection. 

Definition 2.5. (accessibility in finite-dimensional projection) Let C be a 
finite- dimensional sub space of H2(M) andH c be L2-orthogonal projection of 
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i7 2 (M) onto C. The NS/Euler system is time-T accessible in projection on 
C if for any (p G H 2 (M) the image H ( A? ) contains interior points in C. □ 



Now we will introduce notion of solid controllability. 

Definition 2.6. Fix initial condition (p € H2(M) for trajectories of the 
controlled 2D NS/Euler system. Let v(-) € ([0,T];W) be the controlled 
forcing and wt be the corresponding trajectory of the NS system. 

If NS/Euler system is considered on an interval [0,T] (T < +00), then the 
map Et '■ v(-) 1— > wt is called end-point map; the map 11° o T ' /I ~ T is called 
end-point component map, the composition H c o F/T T is called C-projected 
end-point map. □ 

Definition 2.7. Let <E> : M l 1— > M 2 be a continuous map between two metric 
spaces, and S C M. 2 be any subset. We say that $ covers S solidly, if 
S C <I)(A / J 1 ) and this inclusion is stable with respect to C° -small perturbations 
0/$, i.e. for some C° -neighborhood Q of $ and for each map ^ € £1, there 
holds: S C ^(M 1 ). □ 

Definition 2.8. (solid controllability in finite- dimensional projection) The 
2D NS/Euler system is time-T solidly globally controllable in projection on 
finite- dimensional subspace C C H 2 (M), if for any bounded set S in C there 
exists a set of controls Bs such that {l\. c o T '/Tj) (Bs) covers S solidly. □ 

2.4. Problem setting. In this contribution we will discuss the following 
questions. 

• Under what conditions the 2D NS/Euler system is globally control- 
lable in observed component? □ 

• Under what conditions the 2D NS/Euler system is solidly controllable 
in a finite- dimensional projection? □ 

• Under what conditions the 2D NS/Euler system is accessible in a 
finite-dimensional projection? □ 

• Under what conditions the 2D NS/Euler system is L2- approximately 
controllable? □ 

As we explained above the geometry of controllability is encoded in spec- 
tral properties of the Laplacian A and therefore on the geometry of the do- 
main on which the controlled NS system evolves. Below we provide answers 
for particular types of domains. 

3. Geometric control: accessibility and controllability via 



In this section we collect some results of geometric control theory regarding 
accessibility and controllability of finite- dimensional real-analytic control- 
affine systems of the form 



(15) x = f°(x) + f(x)vi(t), x(0) = x°, Vi (t) € R, t = 1, . . . , r. 




Lie brackets 



r 



i=l 
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Geometric approach is coordinate-free, so that it is adapted for dealing with 
dynamics on manifolds, but we will assume that the system (fl~5l) is defined 
on a finite-dimensional linear space M. N in order to maintain parallelism with 
NS systems, which evolve in Hilbert spaces. 

We use standard notation Pt = e** for the flow corresponding to a vector 
field /. 

3.1. Orbits, Lie rank accessibility. Let v(-) S ([0, T]; W) be admis- 
sible controls and x(t) be corresponding trajectories of the system x = 
f°(x) + E[=i P( x )vi(t) with initial point x(0) = x°. We again introduce 
the end-point map £t '■ v(-) i— > x v (T); here x v (-) is the trajectory of ([75]) 
corresponding to the control v(-). 

For each T > time-T (time-< T) attainable set A^ from x° of the sys- 
tem (fT5j) is the image of the set ([0, T];R r ) under the map St, or, equiv- 
alently, the set of points x(T) attained in time T from x° by means of ad- 
missible controls. The time-< T attainable set from x° is ■Ar^ Utg[o XI 
Attainable set from x° of the system (jT5j) is A x o = Ut>o^I"- ^ 

Important notions of geometric control theory are orbits of control system. 

Definition 3.1 (orbits and zero-time orbits of control systems). An orbit of 
the control system / T73]) passing through x° is the set of points obtained from 

x° under the action of ( the group of) diffeomorphisms of the form e* 1 * o- • • o 

e tN f u , where tj G R, j = 1, . . . , N, and f u ' = /°+El=i iS the ri 9 ht ~ 

hand side of / T73]) corresponding to constant control v? = (u\, . . . ,v? r ) G W . 
Zero-time orbit is the subset of the orbit, resulting from the action of these 
diffeomorphisms subject to condition J2j tj = 0. □ 

If we consider 'symmetrization' of the system (fT5l) . 

r 

x = f{x)v +^2f (x) Vi (t), x(0)=x°, v eR, Vi (t) 6 1, i = l,...,r, 
i=l 

then the orbit of (fl~5j) can be interpreted as the attainable set from x° of this 
symmetrization corresponding to application of piecewise-constant controls. 

The famous Nagano theorem relates properties of the orbits and Lie al- 
gebraic properties of the system. It claims that the orbit and the zero-time 
orbit of an analytic system (fT5l) are immersed manifolds of Mr, and tangent 
spaces to these orbits can be calculated via Lie brackets of the vector fields 
{/°,...,/ m }. 

Definition 3.2 (Lie rank and zero-time Lie rank). Take the Lie algebra 
Lie{f°, . . . , f m } generated by {/°, . . . , f m } and evaluate vector fields from 
Lie{f° , . . . , f m } at a point x; the dimension of the resulting linear space 
Lie x {f°, . . . , f m } is Lie rank of the system {/°, . . . , f m } at x. 

Take the Lie ideal generated by spanj/ 1 , . . . , f m } in Lie{f°, . . . , f m } and 
evaluate vector fields from it at x; the dimension of the resulting linear space 
Lie x {f°, . . . , f m }, is zero-time Lie rank at x of the system {f°, . . . , f m }. □ 
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These two Lie ranks either are equal or differ by 1. 

The Nagano theorem claims that in analytic case Lie x {f°, . . . , f m } and 
Lie x {f°, . . . , f m } are tangent spaces at each point x of the orbit and zero- 
time orbit respectively. 

Accessibility properties of analytic control system (fl~5l) are determined by 
its Lie ranks. Recall that the system is called accessible if its attainable set 
A x o has nonempty interior and is strongly accessible if VT > attainable 
sets A^ have nonempty interior. The following fact holds. 

Theorem 3.3 (Jurdjevic-Sussmann (C w case),Krener (C°° case)). If the Lie 
rank of the system of vector fields {f°, . . . , f r } at x° equals n, then VT > 
the interior of the attainable set A~ is nonvoid. Besides A~ possesses 
the interior which is dense in it. If zero-time Lie rank at x° equals n, then 
VT > the interior of the attainable set A^ is nonvoid and this interior is 
dense in A^ - O 

See [20j H for the proof. 

Let £ be a linear subspace of M N and Tl c be a projection of M. N onto C 
A control system (PT5l) is accessible (strongly accessible) from x in projection 
on £ if the image ^\ C 'A x o (IT^^o) contains interior points in £ (for each 
T>0). 

One easily derives from the previous theorem a criterion for accessibility 
in projection. 

Theorem 3.4. IfH c maps Lie x {f°, . . . , f m } (respectively Lie x {f°, . . . , f m } ) 
onto C, then control system {75]) is accessible (respectively strongly accessible) 
at x in projection on C □ 

Proof. The proofs of the two statements are similar; we sketch the proof 
of the first one. Consider the orbit of the system (PT5l) passing through 
xq. The tangent space to the orbit at each of its points x coincides with 
Lie,{/ ,...,/ m }. 

By Theorem 13. 3l the attainable set of the system possesses relative interior 
with respect to the orbit. Moreover there are interior points Xi n t € A x o 
arbitrarily close to x° so that maps Lie xw {/°, . . . , f m } onto C. Then 
sufficiently small neighborhoods of Xi n t in the orbit are contained in A x o and 
are mapped by H c onto a subset of C with nonempty interior. □ 

3.2. Lie extensions and controllability. Controllability is stronger and 
much more delicate property than accessibility. To verify it it does not suffice 
in general to compute the Lie rank which accounts for all the Lie brackets. 
Instead one should select 'good Lie brackets' avoiding 'bad Lie brackets' or 
'obstructions'. 

To have a general idea of what good and bad Lie brackets can be like let 
us consider the following elementary example. 

Example 1. 

X\ = V, X2 = x\. 
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This is 2-dimensional control-affine system (fTSj) with f° = xfd/dx 2 , f 1 = 
d/dx\. The Lie rank of this system equals 2 at each point, the system is 
accessible but uncontrollable from each point x = (x\, x 2 ) given the fact that 
we can not achieve points with x 2 < x 2 . One can prove that the attainable 
set Ax coincides with the half-plane x 2 > x 2 with added point x . One sees 
that that it is possible to move freely (bidirectionally) along good vector field 
/ , while along bad Lie bracket [/ 1 [/ 1 ,/ ]] = 2d/dx 2 we can move only in 
one direction. 

Good Lie brackets form Lie extension of a control system. 

Definition 3.5. The family T' of real analytic vector fields is: 

i) an extension of T if T 1 D T and the closures of the attainable sets 
Ap(x) and A^i{x) coincide; 

ii) time-T extension of T if J 7 ' D T andMT > the closures of the time-T 
attainable sets Ajr(x) and A^p,{x) coincide; 

Hi) fixed-time extension if it is time-T extension VT > 0. 

The vector fields from T \ T are called i) compatible, ii) compatible in 
time T , Hi) compatible in a fixed time with T in the cases i), ii) and Hi) 
correspondingly. □ 

The inclusions Apix) C Af(x), Ajr(x) C ^4^,(x) are obvious. Less 
obvious is the following Proposition (see [2]). 

Proposition 3.6. // an extension T' of an analytic system T is globally 
controllable, then T' also is. □ 

Remark 3.1. When talking about time-T extensions one can consider also 
extensions by time-variant vector fields Xt, t € [0, T] . We will say that vector 
field Xt is time-T compatible with T if it drives the system in time T from 
x to the closure of Ajr(x). □ 

Our idea is to proceed with a series of extensions of a control system 
in order to arrive to extended system for which the controllability can be 
verified and then apply Proposition 13.61 

Obviously the Definition 13.51 is nonconstructive. In what follows we will 
use three particular types of extensions. 

The first natural type is based on possibility of taking topological closure 
of a set of vector fields, maintaining closures of attainable sets. 

Proposition 3.7. (see (20j Ch. 3, §2, Th.5]^ Topological (with respect to 
C°° convergence on compact sets) closure cl{J-) of J- is Lie extension. □ 

The second method underlies the theory of relaxed (or sliding mode) con- 
trols started [El Eli wmc h provided an extension of pioneering contributions 
by L.C.Young [30] and E.J. McShane [23] onto the field of optimal control 
theory. To introduce it we consider a family of so-called relaxation seminorms 
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|| • \\ s k of time-variant vector fields X t , t € [0,T]: 



(16) \\X.\\f K = max 

te[0,T] 



t 

X T \\s,KdT 



with K being a compact in M. N , s > being an integer and ||X,-|| S) _k- being 
a C fc -norm on the compact K. The family of relaxation seminorms define 
relaxation (C°°) topology (metrics) in the set of time- variant vector fields. 

Proposition 3.8 (see |18l 12]). Let a sequence of time-variant vector fields 
X\ converge to X t in the relaxation metrics and let these vector fields have 
compact supports. Then the flows of the vector fields X\ converge to the flow 
of the vector field Xt- □ 

As a corollary of this result one can prove 

Proposition 3.9. For the systems 

( m m \ 

coT = l^Pih, h Gf.fte C"(R N ), /3i > 0, Pi = 1, i = 1,. . . ,m I , 

U=l i=l ) 

and J- the closures of their time-T attainable sets coincide. □ 

Proof of this result and of its modifications can be found in pi Chapter 
8], pOl Chapter 3], [El Chapters II,III ]. 

The next type of extension, we will use, relies upon Lie brackets. It 
appeared in our earlier work on controllability of Euler equation for rigid 
body in [3] and was named there reduction of a control- affine system. Here we 
present a particular version adapted to our problem. Repeated application 
of this extension settles controllability issue for finite- dimensional Galerkin 
approximations of NS systems. 

Proposition 3.10. Consider control- affine analytic system 

(17) ± = f°(x) + f\x)v 1 + f 2 (x)v 2 . 
Let 

(18) [f\f] = 0, [f\[f 1 ,f]]=0. 
Then the system 

x = f°(x) + f\x)n + f(x)v 2 + [f\ [f, /> 12 , 

is fixed-time Lie extension of $17}) . □ 

Sketch of the proof. Take Lipschitzian functions vi(t),V2(t), (fi(0) = 
^2(0) = 0) and change in (fT7l) Vi,v% to e~ 1 i'i(t) + v\ and ev2(t) + V2 respec- 
tively. We arrive to the equation 

(19) x = f°(x) + f\x)v! + f 2 (x)v 2 + (e" 1 / 1 *!^) + ef 2 v 2 (t)) . 
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Applying the 'reduction formula' from [3] or alternatively 'variation of 
constants' formula of chronological calculus ([!]) one can represent the flow 
of (fT9l) as a composition of the flow Pt of the equation 

(20) y = e ad/^)e-V W+ ad/^ 2W/ o (y) + + 
and the flow 

(21) P t = e f 1 e- 1 v 1 (t)+f 2 ev 2 (t) ^ 

For the validity of this decomposition the equality [/ , f 2 ] = is important. 

In (|20j) e a< ^/ is exponential of the operator ad/ : e a< ^/ = X)jLo( a ^/) J ' 
The operator adj is determined by the vector field / and acts on vector fields 
as: ad/<7 = [f,g], [f,g] being the Lie bracket of / and g. 

By first of the relations (TT8l) operators ad/i and ad/2 commute and by the 
second of these relations any iterated Lie bracket of the form (ady^) o ••• o 
(&dp m )f°, with ij = 1,2, vanishes whenever it contains ad/i at least twice. 

Taking the expansion of the operator exponential in (j20j) and using these 
facts we get Q 

(22) y = f°{y) + f\x)v 1 + f\x)v 2 + 
+e- 1 [f 1 ,f](x)v 1 (t) + [f 1 , [fj ]}(x) Vl (t)v 2 (t) + O(£), 

To obtain the flow of (fl~9l) we need to compose the flow of ([22]) with the 
flow (|21~T) . For any fixed T one can get Pt = Id by choosing «].(•)> U2O) such 
that vi(T) = v 2 (T) = 0. 

From now on we deal with the fixed T and the flow of the equation (|22j) . 

Let us change in (122]) Vj(t) to Vj(t) = 2 1 / 2 sm(t/e 2 )vj(t), with Vj(t) being 
functions of bounded variation j = 1,2. The relaxation seminorms of the 
time-variant vector field e -1 ^ 1 , f°](x)2 1 / 2 sm.(t / £ 2 )v\(t) in the right-hand 
side of ((22]) are 0(e) as e -> +0. Also in the right-hand side of (1221) 

[/ 1 ,[/ 2 ,/°]](x)2 S in 2 (t/e 2 )t; 1 (t)^(t) = 

= [Z 1 , [/ 2 ,/°]](a:)«l(t)«2(t) " [A [/ 2 ,/ ]](x)cos(2t/e 2 ^ 1 (t)^(t). 

The relaxation seminorms of the addend [/ , [f 2 , f°]](x) cos(2i / e 2 )vi(t)v2{t) 
are 0(e 2 ) as e — > +0. 

Hence the right-hand sides of the equation (f22l) with the controls = 
2 1 / 2 sin(i/e 2 )'Dj(i), j = 1,2, converge in relaxation metric to the vector field 

f(y) + f\x)vx + f 2 (x)v 2 + [f\ [f 2 , / ]]( S )^(t)^(t) 

as e — > 0. We can consider the product vi(t)v2(t) as a new control V12 and 
invoke the Proposition 13.81 □ 

1 The time-variant vector field abbreviated by O(e) in (f22T) equals 

e0(ead„ 2 (t) f 2 )ad„ 2 (t) /2 f° + e 2 a,d 2 V2 (t) /2 <p{saA V2 (t) p ) [f 1 , f°] , 
where cj>(z) = z^ 1 (e z - 1), ip(z) = z~ 2 (e z - 1 - z). 
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4. Bracket computation in finite and infinite dimensions; 
controlling along "principal axes " 

In this section we adjust the statement of the Proposition 13 .81 for studying 
controllability of the systems fjT2|) and (fl4l) . 

Prom viewpoint of geometric control the Galerkin approximation fj 14[) of 
NS/Euler system is particular type of control-affme system (fl~5l) . Its state 
space is finite- dimensional and is generated by a finite number of eigenfunc- 
tions of the Laplace-Beltrami operator A or modes. 

The dynamics of this control system is determined by quadratic drift vec- 
tor field 

/ o = n o {A-V jg o } + ^A 9 o 

and by controlled forcing f l ( x ) u ii where f l are constant ((^-indepen- 

dent) controlled vector fields. 

Start with computation of the particular Lie brackets, which are involved 
in the formulation of the Proposition 13.81 For two constant vector fields 
f l , f 2 there holds 

Lf, fa] = n° [{f, A-'w} + {w, A- 1 / 4 }) + uAf, i = 1,2; 

[f\ [/ 2 , f o }] = n° ({/ 2 , a- 1 / 1 } + {f\ a- 1 / 2 }) , 
[/ 1 ,[/ 1 ,/ O ]] = 2n°{/ 1 ,A- 1 / 1 }- 

This computation is finite- dimensional but the same holds true if one 
considers constant vector fields acting in infinite- dimensional Hilbert space. 
Taking 'drift' vector field of the equation (|12p in infinite dimensions 

f° = {A- 1 q,q} + uAq. 

we conclude with 

Lemma 4.1. For two constant vector fields f , f 2 there holds 
[f, f] = {f, A-'w} + {w, A- l f} + vAf\ i = 1, 2; 

[/ 1 ,[/ 1 ,/°]] = 2{/ 1 ,A- 1 / 1 }, 

(23) B(f\ f 2 ) = [f\ [f 2 , f ]) = {f 2 , A" 1 / 1 } + {f\ A" 1 / 2 }- □ 

Let us see what is needed for the assumptions of the Proposition 13.101 to 
hold. As long as f , f 2 are constant and hence commuting, all we need is 

(24) [/ 1 ,[/ 1 ,/°]] = {/ 1 ,A- 1 / 1 } = 0. 
In terms of Euler equation for ideal fluid 

^-{A-V^} = 0. 

this means that f 1 corresponds to its steady motion. In particular eigenfunc- 
tions of the Laplace-Beltrami operator A correspond to steady motions and 
satisfy ([24|) . In what follows these eigenfunctions will be used as controlled 
directions. 
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Eigenfunctions of the Laplacian are analogous to principal axes of a (mul- 
tidimensional) rigid body. 

By Proposition 13.101 given two constant controlled vector fields one of 
which corresponds to a steady motion we can extend our control system by 
new controlled vector field ([23]) which is again constant. 

Our method consists of iterating this procedure. The algebraic/geometric 
difficulties, which arise on this way, consist of scrutinizing the newly obtained 
controlled directions. In particular we want to know whether at each step 
we obtain a proper extension. This will be illustrated in the further Sections 
which deal with particular 2D domains. 

Another (analytic) difficulty arises when we pass from finite- dimensional 
approximations to controlled PDE. For these latter the sketched above proof 
of the Proposition 13. 101 is not valid anymore (for example one can not speak 
about flows). We have to reprove statement of the Proposition in each par- 
ticular situation. The main idea will be still based on using fast oscillat- 
ing control and relaxation metric. The analytic difficulties are in proving 
continuity of forcing/trajectory map with respect to such metric. We will 
provide a brief comment in the next Section; the readers can find details in 

[Bid 131 eg. 

5. Controllability and accessibility of Galerkin 
approximations of ns/euler system on t 2 

Here we survey results on accessibility and controllability of Galerkin ap- 
proximations of 2D NS/Euler system on T 2 . 

5.1. Accessibility of Galerkin approximations. The result of the com- 
putation ([231) in the periodic case is easy to visualize when the constant 
controlled vector fields which correspond to the eigenfunctions of Laplacian 
A on T 2 , are written as complex exponentials. 

For two different complex eigenfunctions f 1 = e lk ' x ,f 2 = e ll ' x of the 
Laplacian A on T 2 , where x G M 2 , k,£ E Z 2 , the Poisson bracket in ([23]) 
equals 

(25) B (e ik ' x , e^' x ) = (k A £)(\k\- 2 - \£\- 2 )e^ k+l > x , 

i.e. again corresponds to an eigenfunction of A, provided that \k\ ^ \£\, k A 
£ / 0. The conclusion is that given two pairs of complex exponentials 
e ±ik-x ^ e ±ii-x ag con t r olled vector fields we can extend them by controlled 
vector fields e l ( ±k±i )- x _ 

Iterating the computation of the Lie-Poisson brackets ([23]) and obtaining 
new directions we end up with (finite or infinite) set of functions, which 
contains e ±lk ' x , e and is invariant with respect to the bilinear operation 
S(v). 

Therefore in the case of T 2 starting with controlled vector fields, which 
correspond to the eigenfunctions e lk ' x , k G Z 2 , of the Laplacian, the whole 
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computation of Lie extensions 'can be modeled' on the integer lattice Z 2 of 
"mode indices" k. 

Actually one has to operate with real eigenfunctions of the Laplacian on 
T 2 , i.e. with the functions of the form cos(A; • x),sin(/c • x). Also in this case 
computation of iterated Lie-Poisson brackets j23j) can be modelled on Z 2 
and the addition formulae are similar to those for the complex case. 

We formulate now bracket generating criterion. 

Proposition 5.1 (bracket generating property). // 

(26) + |4 \kA£\ = l, 
then: 

i) invariant with respect to B set of functions, which contains complex 
functions e ±ik ' x , e ±l£ ' x , contains all the eigenfunctions e m ' x , m € 1? \ 0; 

ii) invariant with respect to B set of real functions, which contains 

cos(/c • x), sin(k ■ x), cos(£ ■ x), sin(£ ■ x), 

contains all the eigenfunctions cos(m • x),s'm(m ■ x), m£Z 2 \0. D 

Bracket generating property for Galerkin approximation of 2D NS sys- 
tem under periodic boundary conditions has been established by W.E and 
J.Mathingly in [14]. The following result from [14] is an immediate conse- 
quence of the previous Proposition 15.11 and Theorem 13.31 

Corollary 5.2 (accessibility by means of 4 controls). For any subset M. C 
Z 2 there exists a larger set M' 3 M. such that Galerkin M! -approximation 
controlled by the forcing 

(27) cos(/c ■ x)vk(t) + sin(/c ■ x)wk(t) + cos(£ ■ x)v^{t) + sin(£ • x)wg(t), 

with k, £ satisfying Ii26\) . is strongly accessible. □ 

Here 4 controls Vk(t),Wk(t),V£(t),we(t) are used for providing strong ac- 
cessibility, but actually it can be achieved by a smaller number of controls. 
Example (accessibility by means of 2 controls) Consider forcing 

(28) gv(t) + gv{t), g = cos(k ■ x) + cos{£ ■ x), g = sin(k ■ x) — s'm(£ ■ x). 

Application of controlled forcing (|27|) results in 4 independent controls, each 
one Vk(t),Wk(t),Vi(t),W£(t) entering one and only one of the equations (fl~3l) . 
Under the action of (|28j) each of controls v, v is applied simultaneously to a 
pair of equations (fLS]) . 

Assume ^ \£\, k A £ ^ 0. We compute the bilinear form (|23j) . 

B(g,g) = (-\£\- 2 + \k\- 2 ){cos(k ■ x),cos(£ ■ x)} 

Up to a scalar multiplier B(g, g) equals 

{-\£\' 2 + \k\~ 2 ) sm(k • x) sm{£ ■ x) = 
= (k A £)(-\£\~ 2 + \k\- 2 ) {cos((k -£)-x)- cos((k + £) ■ x)) . 
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Similarly up to a scalar multiplier B(g,g) equals 

(k A £){-\£\~ 2 + \k\- 2 ) (cos((A; - £) ■ x) + cos((A; + £) ■ x)) ; 
The span of B(g,g),B(g,g) coincides with the span of 

g m = cos({k -i)-x), g 21 = cos{(k + £) ■ x). 
The direction 

g m = sin((fc -£)-x) 

is obtained from computation of B(g,g). 

Choose k = (1,1), £ = (1,0); then m = k + £ = (2, 1), n = k - £ = (0, 1). 

Computing new directions B(g 01 ,g),B(g 01 ,g) we note that due to the 
equality \n\ = \£\ they coincide with B(g 01 , sin(/c • x)), i3(^ 01 , cos(/c • x)) re- 
spectively, and their span coincide with the span of functions 

g 12 = sin((fe + n) ■ x)),g 10 = sin((A; — n) ■ x)) = sin(^ • x)) 

Similarly span of the directions 

B(g 01 ,g),B(g°\g) 

coincides with the span of functions 

g 12 = cos((/c + n) ■ x)),g 10 = cos((k — n) ■ x) = cos(£ • x). 

Then g — g 10 = cos(/c ■ x), g — g 10 = sin(k ■ x). 

These two functions together with the functions g 01 ,g 01 form a quadruple, 
which satisfies the assumptions of the Corollary 15.21 Hence our system is 
accessible by means of 2 controls. □ 

Remark 5.1. It is plausible that strong accessibility of Galerkin approxima- 
tion can be achieved by one control. □ 

5.2. Controllability of Galerkin approximations. What regards con- 
trollability, then in general bracket generating property is not sufficient for 
it. One has to select Lie brackets, which lead to Lie extensions, meanwhile 
in the previous example {g, A^ 1 g}, {g, A~ 1 g} a priori do not correspond to 
Lie extensions. 

Even in the finite- dimensional case one needs a stronger result of Proposi- 
tion [3710] to prove controllability property for the Galerkin approximations. 
This has been accomplished in [H [27] in @D and 3D cases. 

Theorem 5.3. Let k,£ satisfy flSj). For any subset M. C Z 2 there exists a 
larger set M.' such that Galerkin M! ' -approximation controlled by the forcing 

cos(k ■ x)vk(t) + sin(A; • x)wk(t) + cos(£ ■ x)vc(t) + sin(£ ■ x)u)£(t), 

is globally controllable. □ 

The proof of this controllability result consists of iterated application of 
Lie extension described in the Proposition 13.101 At each step we extend 
the system by new controlled vector fields which according to fklO) and ([25]) 
correspond to f m±e = cos((m ± £) ■ x), f m±e = sin((m ± £) ■ x). 
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At the end of the iterated procedure we arrive to the system with extended 
set of controls - one for each observed mode. This latter system evidently 
would be controllable. 

An important case when controllability of Galerkin approximation is im- 
plied by bracket generating property regards 2D Euler equation for incom- 
pressible ideal fluid {v = 0), 

Indeed in this case the drift (zero control) dynamics is Hamiltonian and 
it evolves on a compact energy level. By Liouville and Poincare theorems 
Poisson-stable points of this dynamics are dense and one can apply Lobry- 
Bonnard theorem (0[2O]to establish 

Theorem 5.4. Forv = the Galerkin approximation of the 2D Euler system 
controlled is globally controllable by means of the forcing fH[). □ 

We trust that in the case of ideal fluid the controllability of Galerkin 
approximation of 2D Euler system can be achieved by scalar control. 

6. Steady state controlled directions: abstract 
controllability result for ns system 

We can not apply Proposition 13.101 to infinite dimensions directly. Still 
the main idea of adding new controlled directions is valid for NS systems. 
Now we want to formulate an abstract controllability criterion, based on 
Lie extensions and on the computation of Lemma 14.11 This criterion will 
be employed in forthcoming Sections for establishing controllability of NS 
system on various 2D domains. 

Theorem 6.1 (controllability of NS systen via saturation of controls). Let 

span{f 1 ,...,f}=S = V 1 

be a finite dimensional space of controlled directions. Assume f , . . . f r to 
be steady motions of Euler equation (Efy . For each pair of linear subspaces 
C l ,C 2 consider the span of the image B{C l ,C 2 ) of the bilinear map A23\) . 
Define successively 

V j+l = V j + spanB(S j ,V 3 ), j = 1,2,... 

where S J C V 3 is a linear subspace spanned by steady motions. If (J • V J 
is dense in the Sobolev space H 2 (M), then the NS system is controllable in 
finite-dimensional projections and L2- approximately controllable. □ 

As far as V 1 consists of steady motions then T> 3+1 D V 3 + spaDuB(2? , D 3 ). 
Introduce the sequence of spaces 

(29) V{ +1 =V{ + spanB (D 1 , 1^). 

Evidently V\ C V° and density of \Jj V\ in H 2 (M) guarantees controllabil- 
ity. 

Let for f s G V\ D fs = B{f s ,-): 

D fs = {A- l -,f s } + {A- l f s ,-}. 
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The iterated computations ([29]) correspond to iterated applications of the 
operators Df s to f 1 ,. . . f r and taking linear span. 

Corollary 6.2. Let T be the minimal common invariant linear subspace 
of the operators Df x , . . . , Df t which contains fi,...,f k . If J 7 is everywhere 
dense in L 2 (M), then the system is L 2 - approximately controllable and solidly 
controllable in finite dimensional projections. □ 

7. Navier-Stokes and Euler System on T 2 

In this Section we formulate results regarding controllability in finite- 
dimensional projections and L 2 - approximate controllability on T 2 . Namely 
we describe sets of controlled directions which satisfy criterion provided by 
Theorem 16.11 

We take the basis of complex eigenfunctions (e tk ' x ^ of the Laplacian on 
T 2 and take the Fourier expansion of the vorticity w(t,x) = J2 k q k (t)e tk ' x 
and control v(t,x) = J^keic 1 v k(t) e , here k € 71? . As far as w and / are 
real-valued, we have q n = q- n , v n = V- n . We assume Vq = and qo = 0. 

Using the computation (|25l) we write the infinite system of ODE (fl3j) as 

(30) q k = ^2 (m A n)(\m\~ 2 - \n\~ 2 )q m q n - \k\ 2 q k + v k (t). 

m+n=k,\m\<\n\ 

The controls v k are nonvanishing only in the equations for the variables 
q k indexed by symmetric set /C 1 C Z 2 \ {0}. For k K} the dynamics is 

(31) q k = (m A n)(\m\~ 2 - \n\~ 2 )q m q n - \k\ 2 q k , k K} . 

m+n=k,\m\<\n\ 

There is a symmetric set of observed modes K° D /C 1 , which we want 
to steer to some preassigned values. In the only interesting case where /C 1 
is proper subset of JC°, the equations indexed by k S K,° \ /C 1 are of the 
form (|3ip . They do not contain controls and need to be controlled via state 
variables. 

We give a hint of how this can be done; it is an infinite- dimensional version 
of Proposition 13.101 for NS system on T 2 . 

Let r, s G K\ r A s ^ 0, |r| < and k = r + s K . The equations for 
q r ,q s contains controls v r ,v s , while the equation for q k does not. 

Take Lipschitzian functions v r (t),v s (t), (v r (0) = v s (0) = 0) and substitute 
v r , v s in the equations for q r , q s by e~ 1 i' r (t) + v r and ev s (t) + v s . We obtain 
the equations 

q r = (m A n)(\m\~ 2 — \n\~ 2 )q m q n — \r\ 2 q r + e~ 1 v r (t) + v r ; 

m+n=r,\m\<\n\ 

q s = ( m A n)(\m\~ 2 - \n\~ 2 )q m q n - \s\ 2 q s + ev s (t) + v s . 

m+n=s,\m\<\n\ 

Introduce variables q* = q r — e~ 1 v r (t),q* = q s — ev s (t). Assuming v r (T) = 
v s (T) = we conclude q r (T) = q* (T),q s (T) = q* (T) . 
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Let us rewrite the infinite system of ODE (f30]) - ([3"T1) via q*,q* s in place of 
q r ,q s - The right-hand side of the equation for qk = q r + s gains the addend 

(r A s){\r\~ 2 - \s\- 2 )(q; + £~ V(i))(g: + ev s (t)) 

and we see that the controls v r ,v s enter this equation via product v r (t)v s (t). 
Same v r ,v s enter this and all other equations linearly. 

Substitute Vj(t) for j = r,s by Vj(t) = 2 1 / 2 sm(t/e 2 )vj(t) with Vj(t) having 
bounded variations. Then the right-hand side of the equation for qk will gain 
the product 

2sm 2 (t/e 2 )v r (t)v s (t) = (1 - 2 cos(2t/e 2 ))v r (t)v s (t). 

If £ — > +0 this product tends to v r (t)v s (t) in relaxation metric. In all other 
equations v r (t) and v s (t) enter linearly and are multiplied by fast oscillating 
functions. Therefore the corresponding terms tend to in relaxation metric. 

Therefore one can pass (as e — > 0) to a limit system which now contains 
'new' control v rs = v r (t)v s (t) in the equation for qk = q r +s- (This control 
corresponds to the control vu from the Proposition 13 . 10(1 . 

A difficult analytic part is justification of this passage to the limit. It is 
accomplished in [6j [7] for T 2 and in [2U [25] for rectangular and other kinds 
of regular 2D domains. We refer interested readers to these publications. 

Note that the new controlled direction corresponds to complex exponential 
which is an eigenfunction of the Laplacian on T 2 . Hence we can model Lie 
extensions and formulate controllability results in terms of indices k € I? of 
controlled modes. 

Define a sequence of sets KP C Z 2 iteratively as follows: 

(32) j = 2, . . . , 

K? = K?~ x ^J{w- + n\ m,n £ K?~^ j\ ||?n|| ^ \\n\\ A n ^ 0}. 

Definition 7.1. A finite set KL 1 C Z 2 \ {0} of forcing modes is called satu- 
rating if[jf =1 ^ =1?\ {0} ; where are defined by pf. □ 

Theorem 7.2. (controllability in finite- dimensional projection) Let K 1 be a 
saturating set of controlled forcing modes and C be any finite- dimensional 
subspace of H 2 (T 2 ). Then for any T > the NS/Euler system on T 2 is 
time-T solidly controllable in finite-dimensional projections and time-T L2- 
approximately controllable. □ 

As we see the saturating property is crucial for controllability. In [7] the 
following characterization of this property has been established. 

Theorem 7.3. For a symmetric finite set K, 1 = {m 1 , . . . ,m s } C Z 2 the 
following properties are equivalent: 

i) K} is saturating; 

ii) the greatest common divisor of the numbers dij = m l A m? , i,j £ 
{1, . . . , s} equals 1 and there exist m a E K 1 , which are not collinear and 
have different lengths. □ 
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Corollary 7.4. The set K l = {(1, 0), (-1, 0), (1, 1), (-1, -1)} C Z 2 is sat- 
urating. Solid controllability in any finite- dimensional projection and L 2 - 
approximately controllability can be achieved by forcing 4 modes. □ 

8. Controllability of 2D NS system on a rectangular domain 

The study of controllability in finite-dimensional projections and of L 2 - 
approximate controllability on a rectangular domain has been accomplished 
by S.Rodrigues in [24l[25]. The main idea is similar to the one employed for 
the periodic case, but the computations are more intricate. The reason is 
twofold: i) the algebraic properties of the bilinear operation calculated for 
the eigenfunctions of the Laplacian are more complex, and ii) one needs to 
care about boundary conditions. 

For a velocity field u on a rectangular 1Z with sides of lengths a,b (a / b) 
we assume Lions boundary conditions to hold. In terms of the vorticity w 
they can be written as (TTT]) . 

The (vorticity) eigenfunctions (j) k of the Laplacian are 

(33) <j) k = sin Qta) sin (^2^2) , (k = (h,k 2 ) G Z 2 . 

To find extending controlled direction one needs to pick two eigenfunctions 
f 1 = 4> h , f 2 = 4> e , k, t £ Z 2 and to proceed with the computation d23|). This 
results in a linear combination of at most 4 eigenfunctions W s . 

Then again one can follow Lie extensions on two-dimensional lattice Z 2 
of Fourier exponents k = (k%, k 2 ). If the controlled modes are indexed by 
k € K} = {{k%, k 2 )\ 1 < ki, ki < 3, k ^ (3, 3)}, then one can verify that after 
m Lie extensions the set of extended controlled directions will contain all the 
modes (ki, k 2 ) with ki, k 2 < m + 3 with the exception of (m + 3, m + 3). 

This leads to the following controllability result. 

Theorem 8.1 (controllability on rectangular domain). Let 8 controlled di- 
rections correspond to the functions (33]) with k € {(^1,^2)! 1 < ki,k 2 < 
3, k / (3, 3)}. Then the NS system defined on the rectangular domain un- 
der Lions boundary condition is controllable in finite- dimensional projections 
and L 2 - approximately controllable. □ 

Remark 8.1. In the previous Theorem the controllability in projections and 
approximate controllability are guaranteed by 8 controlled directions. □ 

9. Controllability in finite-dimensional projection of 2D 

NS/EULER SYSTEM ON A SMOOTHENED RECTANGULAR DOMAIN 

We are going to prove that the established above controllability of NS sys- 
tem on a rectangular domain persists if one slightly perturbs the rectangular 
and in particular smoothens it. 

To simplify the presentation we assume the rectangular 1Z to have dimen- 
sions 0,6 with a/b irrational. The irrationality of a/b guarantees that the 
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eigenvalues of the Laplacian on 7Z are simple. Let A(0;0), B(0;b), C{a;b), 
D(Q; a) be the vertices of the rectangular. 
Consider parameterized quadruples of points 



Construct e-circumferences centered at these points and round off the cor- 
ners of 1Z by the quarters of these circumferences. We obtain smoothened 
rectangular 7Z £ . 

The following results on continuous dependence of eigenvalues and eigen- 
functions of the Laplacian on the domain holds (see e.g. [12j Vol. I,Ch.VI]): 
for each N choosing e > sufficiently small one can maintain first N eigen- 
values X £ j of A on 1Z £ simple; besides the eigenspace which corresponds to 
j-th eigenvalue A^ of A on 7Z £ converges to the respective eigenspace of A 
on 1Z as well as A^ — > Xj, as e — > 0. 

Some more technical problems arise from the fact that due to boundary 
conditions the evolution space of NS system depends on the boundary and 
one has to take this dependence into account. 

Theorem 9.1 (solid controllability on smoothened rectangular domain). 
Consider NS system on rectangular domain 1Z = [0, a] x [0,6] (b/a - irra- 
tional) controlled along directions which correspond to the eigenfunctions 



Let 1Z £ be the described above perturbations of the rectangle 1Z. Choose for 
each e > eigenfunctions of the Laplacian on 1Z £ ( and respective controlled 
directions which are perturbations of the functions {3$ . Then for any finite- 
dimensional space C there exists ec > such that Ve £ (0, ec] the NS system 
on the perturbed rectangle 7Z £ controlled by these directions is controllable in 
projection on C. □ 

Remark 9.1. Further on there is a need to deal with domains with analytic 
boundary. It is evident that one can slightly perturb the smoothened rectangles 
1Z £ to obtain domains with analytic boundaries (also denoted by 1Z £ ) for which 
the claim of the previous Theorem still holds. □ 

10. Controllability on a generic 2D domain homeomorphic to 



In this section we consider NS system on a 2D domain M under bound- 
ary conditions (ITT|) . We manage to prove that for a generic domain (exact 
meaning of genericity will be specified in a moment) one can choose 8 con- 
trolled directions (corresponding to eigenfunctions of the Laplacian A on 
M) which guarantee controllability in finite- dimensional projections and Li- 
approximate controllability. 

In what follows we assume M to have analytic boundary, to be homeomor- 
phic to a disc and to be endowed with Euclidean metric. Given another such 
domain D one can map M on D analytically (Riemann mapping theorem). 



A £ (s; e),B(e; b — e),C(a — e; b — e 



),D(e;a- e). 



(34) 
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The induced Riemannian metric on D can be represented (see [I3l Vol. 1, 
§§11-13] for an elementary account) in conformal form as: 

(35) v = e a(xi ' X2) (dx 2 1 +dx 2 2 ). 

Its curvature computed as K = (l/2)e~ a ( Xl ' X2 \d 2 a/dx 2 + d 2 a/dx\) must 
vanish, i.e. 

(36) d 2 a/dx\ + d 2 a/dx 2 2 = 0. 

On the contrary if D possesses Riemannian metric \i of form ([35l) . which 
satisfies ([36]) . then D can be isometrically and analytically mapped onto a 
2D domain M' with Euclidean metric. 

Therefore instead of speaking of various 2D domains we will speak about 
various metrics ([35]) satisfying ([36]) on a given domain. 

From now on D will be chosen from set of analytically smoothened rectan- 
gulars for which controllability has been established in the previous section. 
Let TZo = 1Z E0 ; assume that TZo is endowed with a metric ([35]) , for which ([36]) 
holds. Generic domain in M 2 corresponds to a generic metric ([35l) on IZq or 
equivalently to a generic function a(x\,X2) satisfying ([36]) . 

Recall that for the metric ([351) the Laplacian Au is calculated as 

(37) Au = e-< xl ' X2 \d 2 u/dxl + d 2 u/dx 2 2 ). 

Assume that fi,...,fi are eigenfunctions of this Laplacian: A _1 / s = 
\ fs, s = 1, . . . ,1. 

Abstract controllability criterion formulated as Corollary 16.21 claims that 
it suffices to verify that the iterated applications of Df s = {A _1 -,/ s } + 
{A~ 1 f s , •} to fj result in a dense subset of H 2 {M). 

Theorem 10.1. For a residual subset of Riemannian discs M with analytic 
boundaries there exist 8 eigenfunctions (modes) fi, ■ ■ ■ , fs of the Laplace- 
Beltrami operator A on M such that the NS system on M is controllable in 
finite- dimensional projections and L 2 - approximately controllable by means of 
(controlled) forcing applied to the modes fx, . . . , fs. □ 

We call a subset of topological space T residual if it contains an intersec- 
tion of countable family of open dense subsets of T. 
Sketch of the proof: 

For establishing the two types of controllability it suffices to establish 
controllability in projection on any finite- dimensional coordinate subspace 
£ (see [H [6]. Controllability in projection on fixed coordinate subspace £ 
amounts to verification of nonnullity of some determinant, denoted Det£, 
which is calculated via eigenfunctions fi, ■ ■ ■ , fs of A and via their iterated 
Poisson brackets (Corollary 16.21) . 

For the domains 1Z £ constructed in the previous Section and endowed with 
Euclidean metric this determinant Det£ is nonvanishing, provided that e > 
is sufficiently small (say e € (0, eq]). 



24 



A.A.Agrachev, A.V.Sarychev 



Taking any domain M with analytic boundary and Euclidean metric we 
establish its isometry to 1Z £0 provided with metric fj,± = e a ( Xl ' X2 \dxf + dx 2 ) 
which satisfies f[36|) . 

Consider a 'homotopy' 

Ht = e ta{ - xl ' X2 \dx\ + dx 2 2 ), < t < 1. 

All the metrics fit are locally euclidean, given that d 2 a/dx\ + d 2 ajdx\ = 0. 
Obviously the dependence on t of the corresponding Laplacian 

A(t) = e- t < xl ' X2 \d 2 /dxl + d 2 /dx 2 2 ) 

is analytic. 

We want to trace the evolution of a finite number of its eigenvalues A*- 
(j € J- finite set) and of the corresponding eigenfunctions with t varying in 
[0,1]. 

By classical result of perturbation theory (see [211 Ch. 2, Ch. 7]) eigenval- 
ues A* of an analytic family of linear operators are analytic with respect to t 
beyond finite number of exceptional points in [0, 1]. Any moment t for which 
all the eigenvalues A*, (j E J) are distinct is nonexceptional. For generic 
operator singularities of the function t t— » A* may occur when A*- become 
multiple. The eigenvectors and respective eigenprojections may have poles 
at the exceptional points. 

The picture is much more regular for normal operators and in particular 
for Laplacians, which are self-adjoint. In this case the eigenvalues and eigen- 
functions are known to depend analytically on t everywhere on [0, 1] ( |21l 
Ch.2, Th.1.10]). Also the dependence of the derivatives of the eigenfunctions 
on t S [0, 1] is analytic. Hence the determinant Det£ is analytic function of 
t. As far as it is nonvanishing for t = 0, it may vanish only at finite number 
of points t G [0, 1]. 

Take pit corresponding to all nonexceptional t £ [0, 1] for which Det£ is 
nonvanishing. Among those there exist t s arbitrarily close to 1 and metrics 
fit a arbitrarily close to fi\ for which the eigenvalues of interest are distinct 
and Det£ is nonvanishing. The dependence of the eigenfunctions and of their 
derivatives on metric \i is then continuous in a neighborhood of /ij s . Hence 
Det£ is nonzero for all \i from small neighborhoods of these fit s - Taking 
union of these neighborhoods we get an open set whose closure contains fit. 
Repeating the homotopy argument for each [i\ running over the set of all 
metrics l35l)- (|36]) on TZ £0 we get an open dense set of metrics for which Det£ 
is nonvanishing. 

Finally running over the (countable family of all) coordinate subspaces 
C in H 2 {M) and taking the intersection of the corresponding open sets of 
metrics we obtain a residual set of domains for which all the determinants 
DL are nonvanishing and hence the assumption of Corollary 16.21 is verified. 

Remark 10.1. The construction of a residual set of Riemannian metrics 
can be transferred (almost) without alterations to torus T 2 , for which we 
studied controllability of NS/Euler system in the Section^ The conclusion 
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would be that there exists a residual set of Riemannian metrics on T 2 such 
that the assumptions of the Corollary \6.S\ are verified and therefore the NS 
system is controllable in finite- dimensional projections and L 2 - approximately 
controllable on T 2 endowed with any of these metrics. □ 

11. NS/EULER SYSTEM ON SPHERE S 2 

The controlled vector fields we employ in the case of § 2 correspond to 
the eigenfunctions of the corresponding spherical Laplacian or to so called 
spherical harmonics. We start with their brief description. 

11.1. Spherical Harmonics. In this subsection we introduce some notions 
and results regarding spherical harmonics; our source was mainly the book 
[91 Ch. 10,11] by V.I.Arnold. 

Consider sphere S 2 equipped with the Riemannian metrics inherited from 
IR 3 and with area 2-form a. the latter defines symplectic structure on S 2 . 

The eigenfunctions of the spherical Laplacian are described by the fol- 
lowing classical result. Recall that a function g is homogeneous of degree s 
on R n \ 0, if g(nx) = n s g(x) for each k > 0. A function g is harmonic in 
W 1 \ if Ag = 0; where A is the euclidean Laplacian. It is known that a 
harmonic homogeneous function of degree s > is extendable by continuity 
(<?(0) = 0) to a harmonic function on R™. This harmonic function is smooth 
and therefore it must be homogeneous polynomial of integer degree s > 0. 

Theorem 11.1 (|9j). Constants are eigenfunctions of spherical Laplacian 
(of degree 0). If a (smooth) harmonic function defined on M. n \ is homoge- 
neous of degree s > 0, then its restriction onto sphere is eigenfunction of the 
spherical Laplacian A with the eigenvalue —s(s + n — 2). Vice versa every 
eigenfunction of A is a restriction onto S n of a homogeneous polynomial. □ 

Another famous result is the Maxwell's theorem (|9J), which holds in M 3 . 
It states that if p(x) = {x\ + x\ + x 2 ) _1//2 is a fundamental solution of 
the Laplace equation in M 3 , then any spherical harmonic a on S 2 can be 
represented as iterated directional derivative of p: 

a = l± o • • • o l n p, 

where l\, . . . , l n € R 3 and the set {li, . . . , l n } is uniquely determined by a. 

Our controlled directions will correspond to spherical harmonics on S 2 , 
which are the restrictions to § 2 of homogeneous functions on M 3 . In par- 
ticular we invoke so called zonal spherical harmonics, which are iterated 
directional derivatives of p with respect to a fixed direction /. 

Let a, b be smooth (not necessarily homogeneous) functions on M 3 ; the 
Poisson bracket of their restrictions to S 2 can be computed as follows: 

(38) {a\ S 2,b\ s2 }(x) = (x,V x a,V x b), 

where (x, rj, () stays for "mixed product" in R 3 , calculated as the determinant 
of the 3 x 3-matrix whose columns are x, rj, From now on we omit the sign 
of restriction |§2 while writing Poisson bracket. 
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Linear functions (l,x) are, of course, spherical harmonics. We denote by 
I the Hamiltonian field on § 2 associated to the Hamiltonian (l,x), x G § 2 . 
Obviously, I generates rotation of the sphere around the axis I. According 
to the aforesaid la = (x,l,Va) is the Poisson bracket of the functions (l,x) 
and a restricted to S 2 . 

The group of rotations acts (by the change of variables) on the space of 
homogeneous harmonic polynomials of fixed degree n. It is well-known that 
this action is irreducible for any n (see [9] for a sketch of the proof). In other 
words, the following result holds. 

Proposition 11.2. Given a nonzero degree n homogeneous harmonic poly- 
nomial a, the space 

span{l\ o • • • o l k a : k > 0} 
coincides with the space of all degree n homogeneous harmonic polynomials. 

11.2. Poisson brackets of spherical harmonics and controllability. 

Calculating Lie extensions according to the formula ([23]) we obtain iterated 
Poisson brackets of spherical harmonic polynomials, which in general need 
not to be harmonic. 

The following Lemma shows that there is a way of finding some harmonic 
polynomials among them. 

Lemma 11.3. For each n > 2 there exist a harmonic homogeneous poly- 
nomial q of degree 2, and harmonic homogeneous polynomial p of degree 
n > 2 such that their Poisson bracket is again harmonic (and homogeneous 
of degree n + 1) polynomial. 

Proof. Let us take so called quadratic zonal harmonic function q = ^y- 

Being restricted to the sphere S 2 this function coincides with Legendre poly- 
nomial q(x 3 ) = 3x 2 — 1. 

Consider homogeneous harmonic polynomials of variables x±,X2- In polar 
coordinates they are known to have representation r m cos rrup or alterna- 
tively as Re{x\ + 1x2)™, m = 1, 2, We pick the nth degree polynomial 

p(xi, X2) = Re(xi + ix2) n . 

According to (|38l) the Poisson bracket of q,p equals 



{q,p} = (x,S/q,Vp n ) 



xi p' Xl 
x 2 p' X2 
X3 6x3 



-6x3 



x x p' Xl 
x 2 p' X2 
x 3 1 



By (|38j) the latter determinant coincides with e 3 p{x\,X2), where e 3 = 
(0,0,1) is the standard basis vector of M 3 . Hence by Proposition 111.21 the 
value of this determinant is a harmonic polynomial of degree n; it equals 
P(xi, X2) = — x\p' X2 + X2P X1 and therefore does not depend on X3. 

Then {q,p} = — 6x3p(xi, X2). Since both —6x3 and p are harmonic, we 
get A{q,p} = 2V(-6x 3 ) • Vp = -12dp/dx 3 = 0. □ 



Solid Controllability in Fluid Dynamics 



27 



Theorem 11.4. Consider NS/Euler system on sphere S 2 . Let (constant) 
controlled vector fields correspond to three independent linear spherical har- 
monics I 1 , 1 2 , 1 3 , one quadratic harmonic q and one cubic harmonic c. Then 
this set of controlled vector fields is saturating and the NS/Euler system is 
controllable in finite-dimensional projections and Li- approximately control- 
lable. □ 

Proof. It suffices to verify the assumption of the Corollary 16.21 Without lack 
of generality we may think that q = q - the second degree zonal harmonic 
from the previous lemma. Indeed otherwise we may transform q into q by 
taking iterated Poisson brackets with the linear harmonics I 1 , 1 2 , 1 3 . 

In fact taking iterated Poisson brackets of q and c respectively, with 
l\l 2 ,l 3 we obtain all quadratic and cubic harmonics. Thus we manage to 
obtain all the harmonics of degrees < 3. 

Let us proceed by induction with respect to the degree of harmonics. 
Assume that all harmonics of degrees < n are already obtained by taking 
iterated Poisson brackets of {I , I 2 , 1 3 , q, s}. Pick the harmonic polynomial 
p constructed in Lemma 111.31 its Poisson bracket with q is homogeneous 
harmonic polynomial p of degree n + 1, Taking iterated Poisson brackets of 
p with I , I 2 , 1 3 we obtain all polynomials of degree n + 1. □ 

Remark 11.1. Applying the argument similar to the one involved in the pe- 
vious Section one can conclude that there exists a residual set of Riemannian 
metrics on § 2 such that the assumptions of the Corollary \6.S\ are verified and 
therefore the NS system is controllable in finite- dimensional projections and 
L2- approximately controllable on S 2 endowed with any of these metrics. □ 
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